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1,0 INTRODUCTION 


This document presents a solar activity prediction technique for monthly mean 
sunspot numbers over a period of approximately 10 years (from February 1979 to 
January 1989). This includes the predicted maximum epoch of solar cycle 21, ap- 
proximately January I960 and the predicted minimum epoch of solar cycle 22, ap- 
proximately March 198?. Additionally, the solar radio flux 10.7“Centimeter 
smooth values are included for the same time frame using the smooth 13-month 
(sunspot number/ flux) empirical relationship from reference 1. 

The incentive for predicting solar activity values is the requirement of solar 
flux (10.7-centimeter) data as input to upper atmosphere density models utilized 
in mission planning satellite orbital lifetime studies. 

The mathematical analysis used for this technique is discussed in section 3*0. 


2.0 SOLAR ACTIVITY HISTORY AND RELATIONSHIP 

The solar activity sunspot numbers that have been observed and recorded since 
the year 17^9 were obtained from reference 2 and are presented in figure 1 . 

This solar record constitutes a total of 20 solar cycles plus part of the 21st 
cycle, which started in June 1976. These cycles seem to have an approximately 
11 -year cycle that can be identified with the 22-year magnetic period (fig. 2). 

The Zurich listing of data (sunspot number) extends back to the year 17^9; how- 
ever, the reliable data start with solar cycle 8 or approximately 1834, as 
discussed in reference 3 and shown in figure 3. The first seven cycles are 
considered less reliable and were derived from historical data only after the 
Zurich sunspot number was defined. A description of sunspot number is given in 
appendix A. 

Measurements of the 10.7-centimeter solar radio flux have been made since 1947. 
Because of this small data sample, a direct method of predicting the solar radio 
flux has not been developed (ref. 1). A description of solar radio flux is 
presented in appendix B. 
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3.0 M ATHEMATICAL ANALysit 

A mathematical technique employing in part the work of A. K. Paul (ref. 4) 
is given, which allows the reduction of monthly mean sunspot data to a residual 
curve that can bo extropolated, thus allowing a short term prediction. 

Suppose a, b, and h are numbers. Define, for all numbers t, 

F(t) s a cos 2Tiht + b sin Ptrht 

and, as in references 4 and 5, for each number pair f, gj define, for all 
numbers t , 

P(t;f,g) = (exp(2trift) - exp(2Trigt))/iTrt 
Observe that 

P(t;f,g) = 2 sin (iT(f-g)t) exp(27ri(f+g)t/2)/irt 
and also that 


lim P(t;f,g) = 2(f - g) 
t-K) 


Let T(F;f,g) denote 


f P(t;f,g)F(t)dt 


In reference 4, it is shown that 


( 0 


if 

h 

< g 

l(a + 

ib)/2 

if 
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= h 

T(F;f,g) =/ a + 

ib 

if 
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< h < f 

l(a + 

ib)/2 

if 
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< h 
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Let 0 » (a2 + b2)1/2 and notice then tihat there ia a unique number 0 such 
that 

a = 0 cos 0 
and b s c sin 0 

from which it follows that 

F(t) s 0 cos (2Tfht - 0) for all numbers t 

For each number s, define 

Fg(t) s F(t + s) 

= a cos (2TTh(t + 3)) + b sin (2TTh(t+s)) 

= a' cos 2TTht + b' sin 2TTht 

and notice that e = ((a*)^ + As before, there is a unique number 

0 ' such that 

a' = c cos 0' 
and 

b ’ = c sin 0 ’ 

It then follows that 

Fg (t) = c cos (2TTht - G’) 

Recall that Fg(t) = P(t + s); therefore, Fg(t) = c cos(2Trht - (0 - 2irhs)) 
and hence 

0' = 0 - 2irhs 
thus, h = (0 - 0' )/2rfs 
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As before 


T(fs! f,g) r 


( 0 



if 

h < 
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\(a' 


ib')/2 
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An approximation to T(F; f, g) is 


n *~'-oo 


(2/TTDn) sin (ii(f - g)nD) exp(2TTi(f+g)nD/2) F(nD)D 


where D is a positive number. The term for n = 0 is 2(f - g)F(0)D. 

As only a finite amount of data are usually available, we further approximate 
T(F; f, g) by 

N 

(2/it) S (1/n) sin (ir(f - g)nD)exp(2Tii(f + g)nD/2F(nD) 
n=1 

+ 2(f - g) F(0)D 

-N 

+ (2/TJ) (1/n) sin (TT(f - g)nD)exp(2iTi(f + g)nD/2)F(nD) 

n='^1 

Let A = f - g, B = (f + g)/2, D = 1 

F(n) = X„ and for each number s let Tfj(X, s, B, A) be the complex number 
N 

(2/TT) ^ (1/n) sin (TtAn) exp(2TTiBn) Xn^,g 

rl^l 

+ 2(f - g)Xg 

-N 

+ (2/ir) ^ (1/n) sin (TrAn)exp(2TriBn)Xn^.g 

t^-1 
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We now introduce several definitions that ws shall need in order to predict 
sunspot numbers. 

Definition 1 . Suppose c > 0 and f is function defined on 

(,a,bj = (xj a < X < b). The statement that f is c-symmetris means 


|f(x) - f(a + b - x) I < e 
for all X in (a, b). 

Definition 2 . Suppose that R is a function defined on fa, b). The statement 
that R* is the symmetric extension of R to (a, 2b - a) means 


(R(x) if a < x < b 
R“(x) =1 

(R(2b - x) if b < X < 2b - a 

Definition 3 . Suppose v > 0 and F is a function defined on (u, v). The 
statement that F can be e-reduoed on (u, v) means there is a finite sequence 

n 

of periodic functions {gj[} and a finite sequence of nonoverlapping 
m 

intervals {li} i_i filling up (u, v) such that the function R, 
n 

B = F - 81 


is e-symmetric on for i = 1,2...,m. Denote such a reduction by 


I n m ' 

,e,F, (u,v), {gi}i=i , {lihsi, r! 


Notice that there is no claim as to uniqueness of such a reduction. 


Definition 4 . Suppose e > 0 and f is a function defined on (a, b). The 
statement that f is short term e-predictable means 

a. There is an Interval (u, v) containing (a, b) such that b < v, a function 
F defined on (u, v) such that F can be e -reduced on (u, v) and 
f(x) s F(x) for X in (a, b). 


1 

I m 
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b. There ia an e -reduo tion 


n ra 

{e, F, (u, v3, {gi}i=i » {lilisi . R} 


and an Integer j suoh that b is in the right-hand interior of Ij, 

Theorem ; If e > 0 and f is short-term e-prediotable, then f oan be extended 
to a funotion f* suoh that |f*(X) - F(X)| < e where F is the funotion 
stated in definition 

Proof ; Suppose e > 0, f is a funotion defined on (a,b) and F is a funotion 
defined on (u, v) such that F satisfies ooriditions of definition 4. Let 

n m 

{e, F, (u, v), {gilisi, {lilisl; R} 


be a reduction of F satisfying definition k. Then b is in the right-hand 
interior of Ij for_some integer j. Let Ij = (aj, bj) and o the center point 
of (aj, bj). Let R be the restriction of R to (aj, o) and let R* be 
symmetric extension of R* to [aj, 2c - aj) and notice that 2o - aj = bj. Define 


f* s R» + gi on (aj, bj) 

and f» = f 


n 



5 R + gi 

on (a, 

aj3 

Then 

|f»(X) - F(X)j< e 

on (a, 


Suppose 

that we have data 

M 

{Yrlr=0 

at equally spaced time intervals D (where 

D=1), 





0 


2N+1 M 
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Lefc 2N + 1 be an Inbegep less than M and for -N < n < (M - (2N hk 1>) 
let = Xfi+ti 


-N 0 N M - (N + 1) 


Using the notions presented thus far, we have data from (-N, M - (N + 1)) 
and which to predict values for data indexed greater than M - (2N + 1). 

We shall asaitine that there is a function F defined on an Interval (u, v) 
containing (“N, M - <W + 1)3 » which can be e-reduced on (u, v) such that 
for some such reduction definition 4 is satisfied where the f in the definition 
is the known data function. 

n 

We shall generate the sequence of periodic functions {gi}i =1 and the sequence 

m 

nonoverlapping intervals (lillsl by applying the operator T^CX, s, B, A) 
to the data X, 

As we are seeking short-term prediction, we are seeking relatively high 
frequencies between 0.0012 and 0.0833, which correspond to one cycle 
per 5 years and one cycle per month. 

Having the data X, we have three parameters, s, A and B that we can vary. 

The first, s, is a shift, A is a "bandwidth” and B is approximate 
"location." 

Notice that 0 1 a < M - (2N + 1). Let s = 0, select B and a A, then 
compute Tj^(X, 0, B, A), and iterate by reducing A. Then select new B and 
repeat. Continue this process. 

Now select 0 < s < M - (2N + 1) and repeat above process. 

Because D = 1 (1 month in case of sunspot data), we cannot determine frequencies 
that are higher that one cycle per month. Hence 0 < B < 0.0833. 

We continue in this manner, generating six columns of numbers the first being 
always 0, 

3 B A a b c = (a^ + b2)1/2 

We now repeat for a value of 0 < s < M - (2N +1). Hence, we have a subset of 
E^, and we now search for two points within this set 


V-] = (0, B-|, A-], a-j, b-]. Cl) 
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0 33 it first component and the second 

I III I 

Vi s ( s ^ 0, 3i , bi, oi) where ci = oi 

and determine 0i and 0i from which a number hi s (0i - 0i)/2ns is computed. 

Let Ri(n) = - Si(n) for -N < n < N + (M - 2N + 1) 

where 

gl(n) = ai 003 (2nhin) + bi sin (2Tthin). 

A check is now made to determine if X is e-reduced on (u, v) and that the 
reduction satisfies definition M. If not, then we define 

Xin = Ri(n) for -N < n < N + (M - 2N -1) 

Apply now 

T(Xi, 0, BA) 

and T(Xi, s, B, A), and repeat above. 

After each such determination of 

gj(n) = aj cos (2TThjn) + bj sin (2iihjn) 


we set 

Rj(n) :: Rj_i(n) - gj(n) for -N < n < N + (M - 2N - 1) 

where 

Ro(n) = Xn for -N < n < N + (M - 2N - 1) 

At each step Rj is examined to determine if Rj is e-reduced on [-N, N) . 
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Suppose now that we have succeeded In finding an e-reductlon such that defini- 
tion 4 is satisfied. Then by a straightforward application of the theorem, we 
can extend, i.e., predict the future behavior. 


4,0 APPLICATION OF METHOD OF SOLUTION 

By applying the mathematical analysis of section 3*0 to the sunspot number 
data in figure 3> a model was obtained that is composed of 24 spectral lines, 
which were determined by numerous iterations. 

a. Residual curves 

The residual curves (R4, r 8, R12, R16, and R20) obtained by subtracting out 
the frequency, amplitude, and phase (table I) demonstrate the changes in 
the residual curves (fig. 5») 

b. Final residual curve - R24 

Figure 6(a) shows the final residual curve (R24) as a function of 
month/years with the last value at January 1979. Based on the technique 
from section 3.0., this curve was Investigated for short-term symmetry as 
shown in figure 6(b). The shaded portions of the curve indicate symmetry 
from a minimum of 5.3 years to a maximum of 7 years. 

In addition, based on previous tests oases, the final slope indic’sd the re- 
sidual could be extrapolated for 7 or more years. And the data (j idual 
values) are such at the end of the curve the slope has reversed iti® direc- 
tion, therefore making it possible to extrapolate. 

The 20th residual (R20) had the same criteria for extrapolation, except for 
the fact that only a 5- or 6-year prediction of sunspot number would have 
resulted. Therefore, the analysis was not terminated here. Instead, four 
more spectral lines were determined that yielded the final predicted curve 
R24. Figure 7 shows the last 49 years of the final R24 residual curve for 
enlargement purposes. 


c. Determining the predicted residual 

Because the summation of the 24 harmonics of the model with the final resid- 
ual curve (R24) in figure 6 yields the original sunspot record, it was con- 
sidered feasible that extrapolation of the R24 curve an additional 7 to 10 
years results in a prediction of sunspot numbers for 7 to 10 years. 

The method used in the extrapolation is explained as follows: first the 

midpoint of the curve was determined using the digital data of the residual 
(R24) before the last value terminated. By symmetric extension, the resid- 
ual was extrapolated as shown in figure 8. Figure 9 presents the residual 
R24 from 1834 through January 1979 and its 10-year extrapolated residual 
through January 1989. 
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Wlien the 2H harmonics are summed with the curve in figure 8, the original 
sunspot numbers (observed month mean values) and the predicted sunspot num- 
bers are obtained as presented in figure 10. 

NOTE ; The curve shows the sunspot numbers through the maximum epoch of solar 

cycle 21 and minimum epoch of solar cycle 22. The minimum epoch is deter- 
mined when the cycle passes through zero and reverses its algebraic sign. 


5.0 COMPARISON OF PREDICTED SUNSPOT NUMBERS 


A study was conducted to determine the accuracy of using this analysis. The 
same sunspot record was used but terminated at the data point (December 1972). 

The residual curve (R22) used for the final predictable residual is shown in fig- 
ure 11 with the shaded areas for symmetry as in the case of the R24 residual. 

This allowed actual observed sunspot numbers to be compared with the predicted 
results (January 1973 to January 1979) for an accuracy check. The predicted 
sunspot numbers of this study is presented in figure 12. The data are shown in 
the customenary absolute value format (all positive) sunspot numbers (observed 
and predicted) as a function of month/years. The agreement between the two 
sets of data was deemed adequate for an accuracy check. 

a. Sunspot number prediction (February 1979 bo January 1989) 

The 10-year sunspot number prediction from this analysis is presented in fig- 
ure 13 and includes monthly observed data. This plot also shows the predicted 
maximum epoch of cycle 21 and minimum epoch of cycle 22. 

b. Predicted smooth sunspot numbers (13-month smooth means) 

For comparison purposes, 13-month smooth values were computed using the data 
in figure 13. The results are shown with the National Oceanic and Atmospheric 
Administration (NOAA) and Marshall Space Flight Center (MSFC) data in figure 
14. The NOAA data obtained from reference 6 are in the form of predicted 
nominal and positive two-sigma (2 ct) values. The MSFC predicted data (ref. 7) 
are shown in same format. In summation of this plot, the JSC predicted values 
agree comfortably when compared to the NOAA and MSFC (nominal and positive 
2a) values. 

NOTE ; The only way bo truly evaluate these predictions is after several months 
of observed data have been recorded. 


6.0 SUNSPOT NUMBER AND SOLAR RADIO FLUX EMPHERICAL RELATIONSHIP 

One method of determining predicted solar flux values is to use the relationship 
between Ottawa 10 .7-centiraeter flux and Zurich smooth 13-month mean sunspot num- 
bers as used by NOAA in reference 1. This relationship is presented in figure 
15 and was used in this analysis. 
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NOTE: The 10.7-centimeter flux data are used by various organizations as input 

for upper atmosphere density models for orbital lifetime studies. Be- 
cause the solar radio flux is also measured daily, it would be preferable 
to use it directly for predicting solar activity; however, this time se- 
ries is too short (32 years) for a direct prediction technique like the 
sunspot number method. 


7.0 CX) HPARISON OF PREDICTED SOLAP H.'tm u FLUX 

The JSC predicted radio flux (fig. ItJ was derived using the 13-month 
smooth sunspot number history in figure 14 and the flux/sunspot number curve in 
figure 15. This flux data are also included with the predicted NOAA and MSFC 
values, both nominal and plus 2 sigma (refs. 8 and 7, respectively). 

Solar flux predicted data from reference 8 are also included. 

The difference between the predicted flux values shown in this plot is 
attributed to different sunspot number predictions. The predicted flux data 
other than the JSC values were obtained and plotted for this report in August 
1979 and therefore do not represent any updated current predictions since thac 
time. The JSC flux data are presented for a time period of approximately 10 
years, starting in February 1979. 


8.0 CONCLUSIONS 


A method of predicting monthly sunspot numbers over a period of approximately 
10 years is presented using the mathematical analysis technique from section 
3.0, with an additional application step as discussed in section 4.0. The 
monthly mean sunspot numbers used in this study were obtained from reference 2 
and are shown in figure 3. The final predicted sunspot numbers from February 1979 
to January 1989 are presented in figure 13. Also included in this document is a 
comparison of the following solar activity data as follows: 

a. Predicted and observed monthly mean sunspot numbers from January 1973 to 
January 1979 

b. Predicted 13-month smooth sunspot numbers for JSC, MSFC, and NOAA 

c. Predicted solar radio flux 10.7-centimeter data for JSC, MSFC, and NOAA 

Prediction of solar activity sunspot numbers will yield prediction of the solar 
radio flux 10.7-oentimeter values, which are used by various organizations 
as input for upper atmosphere density models for orbital lifetime studies. 
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TABLE I.- SUBTRACTION OF FREQUENCE, AMPLITUDE, AND PHASE 
FROM PREVIOUS RESIDUALS 


R1 = SSN(I) - AMP1(C0S(2tt(FREQ1 x I 
R2 = R1(I\ - AMP2(00S(2it(FREQ2 x I ■ 
R3 = R2(I) - AMP3(C0S(2tt(FREQ3 x I ■ 
R3 = R3(I) - AMP4(C0S(2tt(FREQ4 x I ■ 
RN = r(N-1)(D - AMPN(C0S(2tt(FREQN ; 
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Figure 










Figure 3.- Sunspot number solar cycles (1834-1979), monthly mean value. 
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(a) Epsilon symmetry. 



(b) Symmetric extension. 
Figure 4.- Symmetry. 
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(a) Residual R4. 

Figure 5.- Residuals as a function of time (1834-January 1979) 
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Figure 6.- Residual R24 (1834~January 1979). 
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Figure 9.- Residual R24 and extrapolated residual (1834-January 1989) 
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Figure 11.- Residual R22 symmetry 
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Figure 14.- Predicted smooth sunspot numbers (13-month smooth means) as a function of month/years 








Figure 15.- Sunspot numbers/solar flux, 10. 7- centimeter curve 





Figure 16.- Observed and predicted flux^l0.7-centimeiter values as a function of nionth/years 
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APPENDIX A 

DESCRIPTION OF SUNSPOT NUMBER 

The Zurich relative aunspot number ia an index of the activity of the entire via*- 
ible disk of the Sun. It ia determined each day without reference to preceding 
days. Each isolated cluster of sunspots ia termed a sunspot group, and it may 
consist of one or a large mimber of distinct spots whose size can range from 10 
or more square degrees of the solar surface down to the limit of 1/25 square de- 
gree. 

The relative sunspot number is defined as: 


R s K(10G -f S) 


Where G is the number of sunspot groups, S is the total number of distinct 
spots, and K is a scale factor. 
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APPENDIX B 

DESCRIPTION OF SOLAR RADIO FLUX (10.7 cm) 

Daily observations of the 2800 megahertz radio emissions that originate from the 
solar disk and from any active regions are made at the Algonquin Radio Observa- 
tory (ARO) of the National Research Council of Canada. These are a continuation 
of observations that commenced in Ottawa in 194?. Numerical values of flux are 
published by NOAA in the SGD documents in tables with units of 10“22 wm"2h 2"^ 
and refer to a single calibration made near local noon at 1700 universal time. 
All flux data are then corrected to Sun-Earth distance of one astronautical 
unit . 
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